Domain wall fermion calculation of nucleon g A /g v * 
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We present a preliminary domain-wall fermion lattice-QCD calculation of isovector vector and axial charges, 
g v and g A , of the nucleon. Since the lattice renormalizations, Z v and Z A , of the currents are identical with 
DWF, the lattice ratio (<? A /sv) lattlce directly yields the continuum value. Indeed Z v determined from the matrix 
element of the vector current agrees closely with Z A from a non-perturbative renormalization study of quark 
bilinears. We also obtain spin related quantities Aq/g v and 5q/g v . 



The isovector vector and axial charges, g v and 
g A , of the nucleon provide an interesting addi- 
tional test of the domain wall fermion (DWF) 
method in the baryon sector where it has suc- 
ceeded in reproducing the mass difference be- 
tween the positive- and negative-parity ground 
states, N(939) and iV*(1535) §. These charges 
are defined as 

9v = G f I™ 9v(l 2 ) 

from the isovector vector current {n\V~ (x)\p) = 
iu n [-i ll g v {q 2 ) + qx<T\ f ,g M ( ( l 2 )} u P e ~ lpX , 
and g A = G F lim g A (q 2 ) 

q^—>0 

with the axial current (n\A~(x)\p)= 
iun^^g^q 2 ) + q lJ ,g p {q 2 )]u p e~' ipx . 

The values of g v = G F cos 9 C and g A /g v = 
1.2670(35) are well known from neutron (3 de- 
cay. Here G F denotes the Fermi constant and 
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9 C the Cabibbo angle. g v = G F cos9 c follows 
from vector current conservation. In contrast the 
axial current should receive a strong correction 
from quantum chromodynamics (QCD), resulting 
in the deviation of the ratio g A /g v from unity. 

In lattice calculations in general the two rel- 
evant currents get renormalized by the lattice 
cutoff. With conventional fermion schemes this 
renormalization usually makes the calculations 
rather difficult, if not intractable, even for such 
simple quantities like g v and g A . However with 
DWF it is greatly simplified because Z A — Z v [Q, 
so that the evaluation g^f ttlce / g^ 106 directly 
yields the continuum value. 

Phcnomenological models of baryons have not 
been successful in reproducing this ratio: the non- 
relativistic quark model gives a value of 5/3, and 
the MIT bag model 1.07. Lattice calculations 
typically underestimate g A by 20 % Q. All of 
these previous lattice calculations are done with 
(improved) Wilson fermions and consequently 
suffer from Z A =/= Z v and other renormalization 
complications. 

The present numerical calculations use the 
same gauge configurations reported in ref. H , the 
notations of which we follow here. From this 
work we know DWF works well. In particular: 
1) fermion near-zero mode effects are well under- 
stood, 2) small chiral symmetry breaking induced 
by the finite extra dimension is described by a 
single parameter m res in low-energy effective la- 
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Figure 1. Dependence of N (cross) and N* 
(square and diamond) mass on quark mass, m/. 
Blasts at m/=0 are experimental values in lattice 
unit, a~ 1 ~2 GeV. 



grangian, which decreases as or L s increases, 
with the value of J7i res /m str ango = 0.033(3) at 
(3 = 6.0 and L s = 16, and 3) non-perturbative 
renormalization (NPR) works well for the quark 
bilinears 

Positive-parity nucleon states are created (de- 
stroyed) with interpolating operators By = 



eabc{UaCj5db)u c and B\ 



while the negative-parity ones are created with 
By = 75-B+ = e abc {ulCj 5 d b )j 5 u c and B^ = 
J5B2 — e a b c (u^Cdb)u c with appropriate bound- 
ary conditions in time to reduce backward prop- 
agating contamination By gives the ground- 
state nucleon (iV(939)) mass. On the other hand 
B2 seems to give the first excited positive parity 
state for heavier bare quark mass mf. Whether 
it can reproduce the iV'(1440) mass in the chi- 
ral limit is not yet known. By and B^ masses 
agree with each other, and yield the negative- 
parity ground state, AT* (1535). Our quenched 
DWF calculation reproduces very well this large 
mass splitting between iV(939) and JV*(1535) par- 
ity partners (see Figure [j]). Phenomenological 



models like the non-relativistic quark model and 
the MIT bag model have failed here. It should be 
also noted that an earlier quenched lattice calcu- 
lation using Wilson fermions || failed here too, 
though more recent calculations show improve- 
ments ||. 

So DWF calculation of nucleon matrix elements 
seems promising. g A is interesting because it is 
particularly clean with DWF since Z A =Z V . It is 
also interesting to see how well quenched calcula- 
tions work for a well-known example of soft-pion 
behavior, namely the Goldberger-Treiman rela- 
tion: g A /g v ^if 1T g 1TNN /m N ~ 1.31. We know that 
with DWF the ratio f 7 r/'m N is almost constant 
over the range of m/ we are using, and agrees 
well with the experimental value fa]. 

We follow the standard practice [0] for our two- 
and three-point function calculations. The two- 
point function is defined by 

G jv (i) = Tr[(l + 7t )2{TB 1 (a:)Bi(0)>], 

X 

using B\ — e a bc(u£C'f5db)uc for the proton. The 
three-point function for the local vector current 
is G u /{t,t'), 

Tr[(l + 7t )^(rB 1 (x)^ d (x')i?i(0)}], 

and for the local axial current, (i, £'), 
Tr[(l + 7t ) 7i75 J2{TB 1 (x)A^ d (x')B 1 (0))}, 

x' ,x 

averaged over i = x,y, z. We choose a fixed 
t = ^source — isink and t' < t. From their lattice 
estimates 



G"(t,t')-G d (t,t') 



lattice r 

9t G (t) 

with r = V or A, the continuum values 
g r = Z T g\^, 

are obtained. Here we need the non-perturbative 
renormalizations, defined by 



[uTd] xen = Z r [uTd] 0l 



3 



U™=5 and t sink =21 



i f i Hf**** * J i 



Z v =0.763(5) at m=0.02 
_l i I i I i I i I i I i I i L_ 



4 6 8 10 12 14 16 18 20 22 
Time slice 



0.77 



0.76 



0.75 



0.74 




0.00 0.01 0.02 0.03 0.04 0.05 0.06 



m: quark mass 



Figure 2. Dependence of vector renormalization, 
Z v = l/g^ ttice , on t', at m f = 0.02. A good 
plateau is observed. 



which should satisfy Z A = Z v so that 



continuum 



'G u (t,t')-G d A (t,t') 



lattii 



\Gl{t,t')-G*{t,t>)) 

Note g A is described as Au — Ad where A/ (/ = u 
or d) is defined by 

(p, s\h 5 7^f\p,s) = 2s M Ag, 

with s satisfying s-p = and s 2 = —1. From these 
we obtain spin-polarized longitudinal parton dis- 
tribution, Aq — J dx[qf(x) — q± (x)] = Au + Ad. 
Similarly, df is defined by 

(Pi sN/^75/b, s) = ^{s^Pu - SuP^)Sf, 

with — [7 M ,7 t/ ]/2. This gives the tensor 
charge which is related to the transverse parton 
distribution, 5q = J dx[q±(x) —q-f(x)] — Su + Sd. 
We define G%,(t,t') by inserting T? = qjtJij5Q at 
t' and a projection operator (1 + 74)7175, and 

.lattice = G" T (t,t>)+G d T (t,t>) 

9 G N (t) 

is obtained. Here we need Z T , which is scheme- 
and scale-dependent. Note that Au = 8u = 4/3 
and Ad = Sd = — 1/3 in the heavy quark limit. 



Figure 3. Dependence of vector renormalization, 
Z v = l/g^ ttice , on m/. Note the scale. Slight 
linear dependence extrapolates to the value of 
0.759(6) at m/=0. 



The numerical calculations are from 200 con- 
figurations at (3 = 6.0 on a 16 3 x 32 lattice with 
DWF parameters L s = 16 and M 5 = 1.8. We 
set the source at t = 5, sink at 21, and cur- 
rent insertions in between. The vector renor- 
malization, Z v = l/g|f ttICC , is well-behaved, The 
value 0.763(5) at m, = 0.02 (See Figure |) 
agrees well with Z A = 0.7555(3), obtained from 
(A c ° nse ™ d (t)q l5 q(0)) = Z A (A^(i)q l5 q{0)) [§. 
A linear extrapolation gives Z v = 0.759(6) at 
mf — (Figure [|). For the lattice axial charge, 
giattice^ pi a teaus are seen for 10 < t < 16, with 
a fairly strong dependence on to/ (See for exam- 
ple Figure ||). So the charge ratio, g A /g v , aver- 
aged in 10 < t < 16, linearly extrapolates to 
0.62(13) at to/ = (Figure ||) which is about 
a factor of 2 smaller than experiment. However 
a linear fit may not be justified here. There 
is some curvature apparent in Figure |5|, so the 
value of g A / g v in the chiral limit may be even 
lower. The same calculation yields (with linear 
extrapolations to m/ =0) Au/g v — 0.50(12) and 
Ad/g v = -0.14(6). Similarly, 5u/g v = 0.39(11) 
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Figure 4. The lattice axial charge, g 
rrif = 0.02. A good plateau is seen in 10 < t < 16. 



Figure 5. Dependence of g A /g v on m/. 



and 5d/g v = —0.11(4). A preliminary value for 
Z T /Z A is 1.1(1) §]. 

In summary we have explored the isovector 
weak interaction of the nucleon in lattice QCD 
with domain- wall fermions. All the relevant 
three-point functions are well behaved. Z v deter- 
mined from the matrix element of the vector cur- 
rent agrees closely with that from an NPR study 
of quark bilinears Q]. Linear extrapolations to 
m/=0 give 

. gjg v = 0.62(13), 

. Aq/g v = 0.36(14), 

. 5q/g v = 0.31(12). 

The quite low value of g A / g v that we obtained re- 
quires further investigation. In particular, we are 
studying the Ward-Takahashi identity which gov- 
erns g A . If the matrix element of the pseudoscalar 
density does not develop a pole asm/ — > which 
is expected in the Goldberger-Treiman relation, 
the left hand side, and therefore g A , must vanish. 
Further study is also required to check systematic 
errors arising from finite lattice volume, excited 
states (small separation between Source and t s i n k), 



and quenching (zero modes, absent pion cloud, 
etc), especially in the lighter quark mass region. 
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